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ABSTRACT 


The  study  describes  an  analytical  solution  of  the  nonsimilar 
laminar  boundary  layer  with  pressure  gradient,  variable  wall  tempera¬ 
ture  and  continuous  injection. 

The  method  consists  of  transforming  the  partial  differential 
equations  for  momentum  and  enthalpy  and  then  solving  the  transformed 
equations  by  assuming  polynomial  stream  function  and  enthalpy  pro¬ 
files.  Solutions  obtained  show  very  good  agreement  with  exact 
numerical  results. 

The  solutions  are  obtained  for  flows  over  wedges  as  well  as 
at  the  two-dimensional  stagnation  point  and  over  curved  surfaces 
of  a  two-dimensional  body  in  crossflow.  The  results  of  the  study 
show  that  the  boundary  layer  is  very  strongly  effected  by  the 
injection  mass  flow  rate. 


vix 


INTRODUCTION 


i 


The  concepts  of  transpiration  cooling  with  the  accompanying 
high  efficiency  £sr shielding  a  surface  from  severe  thermal  environ¬ 
ment  are  well  known  and  have  been  discussed  at  length  in  the  litera¬ 
ture.  Critical  temperatures  and  thermal,  environments  are  encountered 
in  combustion  chambers,  over  gas  turbine  blades  and  near  the  nose 
of  a  re-entry  vehicle.  Transpiration  cooling  perhaps  finds  its  most 
important  applications  in  keeping  these  components  at  an  acceptable 
temperature. 

The  advantages  of  increased  gas  temperature  on  the  performance 
of  a  gas  turbine  are  well  known.  'The  increase  in  turbine  entry 
temperature  causes  an  almost  monotonic  rise  in  efficiency  and 
specific  power.  On  the  other  hand,  this  high  temperature  coupled 
with  excessive  centrifugal  and  gas  bending  stresses  produce  a 
critical  situation  in  the  gas  turbines,  especially  over  the  first 
stage  rotor  blades.  To  insure  efficient  performance  over  a  reason¬ 
ably  high  lifetime  makes  the  problem  of  blade  cooling  absolutely 
necessary. 

In  a  previous  report ^ ,  the  authors  have  described  an 
analytical  solution  of  the  compressible  laminar  boundary  layer 
over  a  flat  plate  with  continuous  or  local  mass  injection.  There, 
the  authors  showed  that  even  a  small  injection  rate  (2%  or  less 
of  the  main  flow)  is  very  effective  in  lowering  the  plate  wall 
temperature.  Though  high  pressure  gradients  are  present  near  the 
stagnation  point  and  over  the  rest  of  the  turbine  blade,  similar 
results  are  observed  for  this  case  also.  For  example,  Gollnick'2^ 
showed  that  transpiration  cooling  is  sufficient  for  practical  use 
in  cooling  the  axisymmetric  stagnation  region  of  a  re-entry  vehicle. 
Goodwin  and  Howe'3)  showed  that  for  relatively  low  mass  flow  rate 
(0.1%  to  1.0%  of  the  main  stream)  the  aerodynamic  heating  is 
shielded  to  a  great  extent.  These  solutions,  however,  are  based 
on  the  principle  of  "similarity"  and  have  to  be  viewed  critically. 

In  the  presence  of  injection,  the  boundary  layer  is  not  likely 
to  be  similar.  Various  authors (4) ,  (5) ,  (6)  have  solved  these  non¬ 
similar  boundary  layers  using  analytical  or  numerical  techniques. 
These  analytical  solutions  predict  nonsimilar  boundary  layer 
behavior  well  for  uniform  injection  along  a  flat  plate  with  con¬ 
stant  fluid  properties.  However,  these  approaches  breakdown  in 
the  presence  of  pressure  gradients. 

This  study  is  an  extension  of  the  previous  work ^  and  in¬ 
corporates  the  presence  of  pressure  gradients  in  the  laminar  flow 
field.  The  report  describes  an  analytical  solution  of  the  non¬ 
similar  laminar  boundary  layer  with  continuous  injection  and 
variable  wall  temperature  over  a  two-dimensional  body  in  cross- 
flow.  The  method,  which  consists  of  transforming  the  partial 
differentia]  equations  of  momentum  and  enthalpy  and  then  solving 
them  by  assuming  polynomial  velocity  and  enthalpy  profiles,  is 
similar  to  the  one  described  in  the  previous  report'1). 
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The  aerodynamic  performance  of  the  heavily  loaded  turbine 
blade  necessitates  that  the  mass  flow  rate  of  the  coolant  effusing 
from  the  blade  surface  be  kept  as  low  as  possible.  First,  the  low, 
transverse  velocity  of  the,  injected  coolant  will  cause  a  decrease 
in  the  kinetic  energy  of  the  mainstream  through  direct  momentum 
exchange.  Second,  as  shown  in  the  first  report,  the  injection  in 
the  laminar  boundary  layer  makes  the  layer  unstable  and  enhances 
thetransition  to  turbulence.  Consequently,  the  injection  flow  rate 
is  to  be  kept  as  low  as  possible.  In  this  study,  the  injection 
rate  is  assumed  1%  or  less  of  the  main  flow  rate  near  the  stagna¬ 
tion  point  and  less  than  5%  over  the  rest  of  the  body  surface. 


ANALYSIS 


The  equations  of  motion  for  a  two-dimensional  body  in  cross 
flow  are : 

Continuity; 


!_  (pu)  +  ~  (pv)  =  0  (1) 

3x  3y 


Momentum; 


--  3u  --  3u  3p  3 

pU  — —  +  pV  +  ( y 

3x  3y  3x  3y 


3U, 


9y 


(2) 


Energy; 


--  3  H  ,  --  3H 


M  3H  ,  -  .. 
lPj«  _  V*  (i 
3x  3y  3y  3y 


pu  —  +  pv  —  =  — —  [ 


l_,  1  3u2  -I 
Pr'  o'-  i 

3y 


(3) 


Here_it  is  assumed  that  the  flow  is  subsonic, 
enthalpy  H  is  given  by 


The  total 


H  =  C  T  +  i  u2 
P  2 


(4) 


(11) 


and  the  various  thermodynamic  and  transport  variables  reduce  to 


1  2 

T  _H-fu 

Te  H-  -  i  u  2 
3  2  e 


(12) 


Here  a  linear  viscosity- temperature  relationship  is  assumed. 

It  is  convenient  at  this  point  to  introduce  the  Lees-Dorodnitsyn 
transformation  of  coordinates. 
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C  -  /  Pe"eUe  dx 


(13) 


and 


peye  y  o 

n  =  - -  /  ^ —  dy 

✓25  o  pe 


The  stream  function  \p  is  defined  as 


3l|)  _ 
3  X 


PV 


and 


it  - 

3y  “ 


pu 


and  is  assumed  to  have  the  form 


4<  =  m  f  U  ,n ) 


(14) 


(15) 


(16) 


From  Equations  (14),  (15)  and  (16) 

(5,n)  (17) 

e 

where  (  ) 1  denotes  differentiation  with  respect  to  n 
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An  enthalpy  function  is 'defined  as 


H 

ir  =  g 


(18) 


The  transformed  fundamental  equations  obtained  from  the 
above  relations  are  now  written  as 


f.t,  +  ff.«  =  25  [f1  f"  ||]  +  X  (g  -  f'2) 


(19) 


and 


g"  +  fg‘  =  2?[f«  f|  -  g*  ff] 


(20) 


where 


.  .  . 

ue 


(21) 


BOUNDARY  CONDITIONS 


The  boundary  conditions  in  the  physical  plane  are- 


u  =  0 


-  _  o  ^  (pv)w  =  w  =  constant 


E  7?  -  Gv)w  cp  (tw  -  tc) 
°y 


(22) 


u  =  U 


T  =  T 


(23) 


Third  condition  of  Equation  (22)  means  that  all  the  h«_at 
transferred  from  the  hot  gas  goes  to  increase  the  coolant  temp¬ 
erature  to  the  exit  wall  temperature.  In  terms  of  the  trans¬ 
formed  variables,  these  boundary  conditions  reduce  to. 
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n  + 


n.  =  0  { 


'  V  = 

0 

i 

w  ^Re 

f  = 

o 

w 

Peye  ^ 

o  ue 

2fl  = 

/2§”  w  /ReQ 

{^w 

k.  3n 

u 

e 

\  V  - 

1 

g  =  1 


{24) 


(25) 


The  second  condition  of  Equation  (24)  is  derived  from  the 
definition  of  the  stream  functions.  Equations  (15)  and  (16). 

For  flow  over  wedges  where  the  pressure  gradient  (X)  is  constant, 
Ue,can  be  written  as  an  explicit  function  of  C.  In  that  case, 
the  value  of  fw  can  be  obtained  as  an  explicit  function  of  £. 
However,  for  flow  over  curved  surfaces,  pressure  gradient  varies 
with  the  downstream  distance  and  hence  it  is  much  easier  to  obtain 
the  value  of  ^w  at  each  point  (?)  by  numerically  carrying  out  the 
integration  up  to  that  location. 


SOLUTION  OF  THE  NONSIMILAR  EQUATIONS 


The  solution  of  the  nonsimilar  Equations  (19)  and  (20)  is 
obtained  in  a  manner  analogous  to  that  described  in  the  previous 
report^1).  The  modified  stream,  and  the  enthalpy  function  are 
expanded  in  terms  of  polynomials  of  n : 


f 


9 

N  +  E 
n=l 


a  n 
n 

n! 


n 


(26) 


g 


7 

+  E 
n=l 


b  ri 
n 


n 


n! 


(27) 


where  a  and  b  are  functions  of  £  only, 
n  n  -1 

From  previous  work^'^  it  was  realized  that  to  obtain  a 
reasonably  good  accuracy  neai  the  two-dimensional  stagnation  point 
(X  =  -  1.0),  it  is  necessary  to  include  at  least  10  terms  in  the 
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(7) 


stream  function.  It  was  also  observed'"  that  the  momentum  equation 
is  much  more  sensitive  and  unstable  than  the  energy  equation.  Con¬ 
sequently,  terms  up  to  are  included  in  the  momentum  equation 
while  energy  equations  retain  terms  only  up  to  n^. 


If  all  the  terms  present  in  the  coefficients  of  powers  of  n 
are  retained,  the  solution  becomes  very  lengthy  and  cumbersome. 
Hence,  an  order  of  magnitude  analysis  is  necessary  to  retain  only 
the  most  significant  terms. 


From  the  boundary  conditions  at  n  =  0  and  from  direct  dif¬ 
ferentiation  of  the  Equations  (26)  and  (27)  it  follows 


a!  =  0 


N  =  f 


w 


b,  =  fa 
1 


(28) 

(29) 

(30) 


1=0 


Again,  following  the  procedure  described  for  the  flat  plate 
solution,  the  unknown  coefficients  a3 ,  a^j  b2,  b3  etc.  are  obtained 
in  terms  of  a2,  N  and  gw  by  substituting  Equations  (26)  and  (27) 
into  Equations  (19)  and  (20)  and  matching  the  coefficient  of  powers 
of  n* 


Since  the  injection  mass  flow  rate  is  kept  1%  or  less  (of 
the  main  flow),  the  absolute  value  of  N,  especially  near  the 
stagnation  point  turns  out  to  be  always  less  than  0.7.  In  that 
case,  the  contribution  of  all  terms  involving  N  in  a~j ,  as,  ag 
and  bg  and  b7  is  quite  small  as  compared  to  the  contribution  of 
the  terms  present  in  a2  to  ag  and  b]_  to  b5»  Hence  all  terms 
involving  N  are  neglected  in  the  expression  of  a-j ,  as,  ag,  bg 
and  b7.  Again,  only  first  and  second  order  terms  in  N  are  re¬ 
tained  in  ag  and  b^.  Also,  since  the  values  of  ?  da/d?,  ?  dgw/d?, 

?  dN/d?,  b2  etc.  are  the  same  order  of  magnitude  (or  less)  as  N, 
the  same  assumptions  are  applied  to  these  terms  or  a  product  of 
these  terms.  These  assumptions  greatly  simplify  the  final  solution. 
Further  comments  on  these  assumptions  are  made  in  later  chapters. 

The  various  coefficients  are  obtained  as 


=  a  =  unknown 


(31) 


a0  =  -  Na  -  2?  a  (N)  _  +  X  a 


(32) 


„2 

N  a 


45aN (N)  ^  +  4S2a(N)2  -  NAct  +  bxA 


2?Xa(N)? 


>‘%V' 


-  N3a  -  6£aN2  (N)  -  12£2aN(N)2  -  a2  +  2a  (a)  _  -  8£3a(N)3 

4  5  5  5 

+  XN2a  -  bxNX  +  4£AaN(N)^  -  2Xa2  -  2?b1X{N)^  +  b2X 
+  4?2Act(N)  (34) 


ag  =  -  Na5  -  4a3a  +  Xb^  -  6Aaa3  -  2?(N)^a5  +  4?a(a3)^ 


Neglecting  all  terms  of  the  order  of  w  or  higher  and  sub¬ 
stituting  the  values  of  a3,  a5,  b3,  etc. 

a6  =  3N2b1X  +  5Na2  +  8ANa2  -  65aN(a)?  +  12Xb1N5(N)? 

-  4aXa  +  2£aA (a)  +  12?2b.X(N)2  -  4£2a (N) f (a) _ 

5  l  5  5  £ 

+  16Aa2£  (N)  r  -  6X2aa  +  45aa(A)_  +  6agA(a)._ 

5  5  5 

-  2  a2£(N)_  -  852a  (a)  (N)  (3* 

5  5  5 


a7  =  -  Nag  -  7aa4  -  4(a3)2  +  Xb4  -  6X(a3)2  -  8Aaa4 


+  6?a(a4)?  -  25a4(a)?  +  25a3(a3)?  -  25ag(N)^ 


This  reduces  to 


a?  =  -  8ab1X  -  4X2a2  -  6X3o2  -  8X2ab1 


r*; 
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Similarly,  expressions  for  ag  and  ag  are 

a  =  -  16X2bxa  +  32a3X  +  20a3X2  +  11a3  -  20X3b1a 


=  90a2Xa  +  226X2a2a  -  lSb^X2  -  20b12X3  +  132X3a2a  (40) 


Here 


a  =  a 


b  S12, 

°1  30 


b2  =  “  Nb1  -  25b1(N)? 

b3  "  w2bl  +  4?Nbi<N>£  +  2a(a)?£+  45^(1$) 2 
b4  =  -  ab1  -  N3b1  -  6N2b1(N)^-  4?Na(a)? 

-  1252Nb1(N?^  +  458^)  ?  -  8$2a(a)  ^ 

-  25b1(a)5  ~  8?3b1(N)3  +  25Xa(a)? 

b5  =  5ab1N  -  165aN(b1)?;  -  4CNXa(a)^  -  Xab1 
+  lO^b^N  (a)  +  14Cab1(N)?  +  6?Xa(b1)^ 

-  16£2Xa  (a)  ^  (N)  ^  -  48£2a  (b-^)  ^  (N)  ^  +  285  (a)^bi^)^ 

-  2^0  (X)  +  (a)e(65N2a  +  5653a(N)2  +  4052aN(N)c> 


Similarly, 


b6  =  -  bi  X 
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b?  =  llb1a2  +  2bxXa2 


(49) 


All  other  terms  in  b .  and  b~  are  small  in  comparison  to  the  above 
terms  and  are  neglected. 

The  momentum  equation,  after  first  integration  is  written  as 
f' '  =  e"?  w(5,n)  (50) 


where 

n 

P  =  /  f  dn  (51) 

o 


[J 

y 

u 

LI 

u 


and 


u(5,n)  =  /  [2C{f* (f • ) 


o 


f"  (f)  }  +  (g  -  f'2)]  eF  dn  +  a 

(52) 


Now  to  (?, n)  is  expanded  in  terms  of  powers  of  n  as 
oj  n 

u (£, n)  =  2 
o 


d  n" 
n 

ni 


(53) 


Equations  (26)  and  (53)  together  with  the  well  known  expansion 
of  the  exponent 


etc. 


(54) 


are  then  substituted  into  Equation  (50).  Again,  matching  the 
coefficients  of  powers  of  n  gives  the  coefficients  "d  "  in  terms 
of  a,  a  and  N.  n 

A  second  integration  reduces  Equation  (50)  to 

n  „  “  d  nn 

f*  (5,n>  =  f  e  (E  ~y~)  dn  (55) 

o  o 


i 
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As  shown  elsewhere 


U) 


r 


this  equation  can  be  written  as 


n  -  e 

f' (S,n)  =  /  e  6  (  i  I  nn)  dn  (55) 

o  n=0  n 


Here  the  right  hand  side  is  expanded  up  to  the  eighth  power  of  n . 
ln  are  functions  of  cL,  and  hence  functions  of  a,  a  and  N. 

By  setting  the  upper  limit  to  n  =  «  gives  the  left  hand  side 
as 


f'  (5,»)  =  1.0 


(57) 


(This  follows  from  the  remaining  boundary  condition.)  Thus,  when 
the  integration  is  performed  and  the  numerical  values  of  the  Gamma 
functions  substituted,  the  final  solution  of  the  momentum  equation 
is  obtained  as 


a2/3‘Co+  720  C6  -  X<5  C3  C6>  +36  V2)  +  a1/3[N(-C1+  ij  ^ 

'*  T  ic4>  +  «<N>5<-2CX  -  12  C4  +  I  1C4))+[2  C2«2(N)2 
+  2  C2(»(»l£  +  |  C2  N2]+  a‘1/3[§  C3  5(a),  -  C,5(a> 


-3  3 


3  C35j(N)5  -  2  C3C  N(N)^  -  C3  IT  (N)  -  ~  NJJ 


£  360  6  '**'5 

C3  ,,3. 


•2/3 


1  ,2, 


+  a  ■'^[a(XCi“I  X‘C4)  +  Jc45A(a)e  +  «C7(^+A_+t|_  x3) 

-  J  C4N£(a)?-  |  C452(N)?(a)?+  |  }  +  a”1  [a  (0.5  C^X 

-  C25X(H)?  +  ij.  CSNX+  C55X(N)C  +l  C5NX2)  +  bl(|  c2x  C5x 

-  X5  c5»2)-  J  C55NX(a)5  +  |  ?N  C^ol  (X )  J  +  a"4/3  [a  (|  CjXN2 

+  |  C3N5X(N)5+  |  c3c2x(n)2)+  to1  (-  |  c3ha-  |  c3?x(n) 


+  I?4  CSNX  +  15  C6«X  (H>  5  +  15  °6N>2> )  +  a"5/3  Ibl  <|  C4n2x 
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5 

30 


5 

20 


+  |  C4ASN(N)?  +  ~  C4§2X(N)^)]+  a“2E«2(x§60  “  ■ 

,  17  n  .3  .  7  .  4  _  ,  ,  -7/3  r13  ,  .,,2  2  1  „  .  ,2 

+  2520  C8X  +  840  X  C8^  +  ^240  C6NX  °  ”  45  CSblX  ° 

+  125"  a2“2W)?  +  IJq  C6NX3u2]  (58) 


where  C  =  1.623,  Cx  =  1.493,  C2  =  2.0,  C3  =  3.234,  C4  =  3.73, 

C5  =  1270,  Cg  =  25.92,  C?  =!59.6  and  Cg  =  144.0. 

( 8 ) 

As  shewn  by  Meksyn  ,  the  solution  of  the  momentum  equation. 
Equation  (58) ,  may  or  may  not  be  convergent  depending  on  the  number 
of  terms  involved  and  the  value  of  X.  Setting  N  =  (N)r  =  (a)^  =  (a)^  = 
and  a  =  1,  it  is  possible  to  evaluate  the  values  of  "a"  for  various 
values  of  X.  These  values  are  presented  in  Table  I. 


Table  I 

Values  of  Shear  Stress  Parameter  as  a  Function  of  X  (from  Equation  (58)) 


X 

aLt 

Literature 

Values 

dpr 

without  Euler 
Trans  formation 

a  „ 
pr 

with  Euler 

Trans  formation 

-1.0 

1.232 

1.31 

1,234 

in 

• 

o 

l 

0.927 

0.94 

0.933 

i 

o 

• 

to 

0.687 

0.696 

0.694 

1 

o 

• 

l-1 

0.587 

0.592 

0.588 

0 

0.47 

0.48 

0.481 

Comparing  the  values  of  apr  (without  Euler  transformation  from 
Table  I  and  from  previous  work  (D ,  (7)  at  X  -  1.0,  it  is  evident 
that  at  that  point  the  solution  is  very  slowly  converging. 

Meksyn  ;  has  shown  that,  this  is  due  to  the  limited  radius  of 
convergence  of  Equation  (53).  To  improve  the  convergence,  Meksyn 
applies  the  Euler  Transformation. 


It  is  necessary  to  realize  that  the  expression  for  u(5,r, ), 

Equation  (52)  itself ,  is  convergent.  It  is  the  transformed  series. 
Equation  (53)  that  is  very  slowly  converging.  As  shown  by  Euler, 
the  Siam  of  this  very  slowly  converging  series  is  the  finite  numerical 
value  of  the  convergent  expression.  Equation  (52).  A  detailed  ac- 
count  of  the  Euler  transformation  is  given  by  Bromwich  and  Macrobert^  , 
Hardy (10)  and  Knopp(H).  Applying  Euler  transformation  to  the  prese.nt 
solution,  the  values  of  the  coefficients  become  CQ  =  1.621,  =  1.464, 

C2  =  1.82,  C3  =  2.41,  C4  =  2.98,  C5  =  3.05,  Cg  =  2.33,  C7  =  ltl63  and 
Cq  =  0.281.  Substituting  these  values  in  Equation  (58)  and  again 
setting  N  =  (N)r  =  (a)£  =  (a)r  =  0  and  a  =  1,  the  values  of  "a"  as 
as  a  function  of  X  are^re-evaluated.  These  values  are  also  shown  in 
Table  I.  These  values  of  "a"  obtained  with  Euler  transformation 
agree  extremely  well  with  the  literature  values. 

In  any  case,  from  Table  I,  it  can  be  concluded  that  the  Euler 
transformation  gives  significant  improvements  only  for  high  values 
of  X,  i.e.,  near  X  =  -1.0.  For  cases  where  X  lies  between  X  =  0 
and  X  =  -0.7,  the  solution  of  the  momentum  Equation  (58),  is  con¬ 
vergent  and  the  improvement  caused  by  the  Euler  transformation  is 
small. 

As  will  be  seen  later,  final  solution  of  the  momentum  equation 
is  obtained  with,  as  well  as,  without  the  Euler  transformation. 

Comparison  between  these  two  solutions  is  made  in  the  chapter  "Results 
and  Discussion". 

Solution  of  the  energy  equation  is  obtained  along  similar 
lines,  with  the  final  solution  being 

45(a)?  +  a'-^fl.esgbj^  +  0.072b1X  -  6.476Nf(a) 

-  4.312£2(a)5(N)?  +  5.0365(b]_)?] 

+  a"2/3[-  3. 9711  (b,)  +  £(a)  (1.49N2  +  8.945N(N) 

*  s  s  s 

+  13.92c2(N)2)  -  11. 92£2  (b, )  (N)  -  2.9865b,  (N) r  -  1.245Nb1 

+  0.4965b, x  (N)  1  +  a"1 [452bn (N) 2  +  45Nb,  (N)  +  b,N2] 

15  15  I  5  f 

+  a"4/3[2.518X5a(a)c-  2 . 5 89 5bx (a)  ]  +  a“5/3  [-0. 248Xbxa  + 2 . 432 5Nb1 (a) 

-  0.9935NXa(a)  +  1.495Xa(b,)_  -  3. 9752Xa  (a) c  (N)  +  6.9552b1  (N)  (a) 

5  1  5  55  155 

-  0.497cbia(X)  ]  +  a~2[-0.6b1aXN  -  0.1b12X  -  3NXb]a(l  +  X)/2 
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+  N^a  (X)  0.2b1?Xa  (N)^]  +  a~7//3  [0.288b12NX  +  0.144^-^  (N)  ?]-  1-  a 

(59) 


Regardless  cf  the  number  of  terms  involved  or  the  value  of  X, 
Equation  (59)  will  always  give  a  =  1.0  at  the  stagnation  point 
for  N  =  0.  Hence,  Euler  transformation  is  not  applied  to  the 
solution  of  the  energy  equation. 


Conditions  at  the  Stagnation  Point 

To  solve  the  Equations  (58)  and  (59)  simultaneously  it  is 
necessary  to  have  the  initial  values  of  "a"  and  a  i.e.,  a^n  and 
ain*  These  initial  values  are  also  stagnation  point  values.  There 
are  two  sets  of  stagnation  point  values  depending  on  the  mode  of 
injection. 

(i)  If  the  injection  is  assumed  to  occur  from  x  >  0,  the  stag¬ 
nation  values  of  a  are  the  same  as  those  given  in  Table  II,  i.e., 
for  X  =  -  1.0, 


a^n  =  1.234  (with  Euler  transformation) 

a.  =  1.31  (without  Euler  transformation) 


and 


a 


in 


1.0 


(ii)  If  the  injection  is  assumed  to  occur  at  the  stagnation 
point,  the  problem  can  be  resolved  as  follows. 

Very  near  the  stagnation  point,  the  velocity  outside  the 
boundary  layer  of  a  two-dimensional  body  in  crossflow  can  be 
written  as 


U  =  £x  (60) 

e 

where  £  is  a  constant. 

For  a  two-dimensional  elliptic  cylinder  in  crossflow, 


£ 


1  +  A 


(61) 
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where  A  is  the  ratio  of  minor  axis  to  major  axis.  From  the 
definition  of  K,  and  Equation  (60), 


(62) 


Again,  for  very  small  values  of  £da/d£,  gda/d E,  and 
£dN/d£  are  neglected.  Then, 


(63) 


Equation  (63)  then  gives  the  stagnation  value  of  N.  The 
stagnation  point  values  of  "a"  and  a  can  now  be  obtained  from 
simultaneous  solution  of  Equations  (58)  and  (59)  using  Equations 
(60)  and  (63).  As  mentioned  above,  £da/d£,  £da/dS,  £dA/d£  and 
£dN/df;  are  omitted.  Under  these  assumptions,  the  solution  reduces 
to  "similar  solution".  However,  at  the  stagnation  point  the  value 
of  £  ->  0  and  the  nonsimilar  terms  tend  to  drop  out. 

Since  the  factor  fw(N)  incorporates  the*  main  flow  rate,  body 
geometry  and  the  injection  flow  rate,  the  stagnation  point  solutions 
obtained  with  this  factor  as  a  parameter  would  constitute  a  very 
"general"  solution.  That  is,  such  solutions  would  apply  to  any 
two-dimensional  body  in  crossflow.  The  stagnation  point  values  are 
given  in  Table  II. 


Table  II 


A 

"fw 

a 

(1  -  a)/(  1  -  Tc) 

-1.0 

0.0707 

1.16 

0.113 

i 

H* 

• 

o 

0.212 

1.02 

0.299 

-1.0 

0.354 

0.885 

0.456 

-1.0 

0.707 

0.556 

0.783 

Once  the  value  of  "a"  and  of  a  has  been  determined,  the 
friction  factor  and  the  heat  transfer  coefficient  are  obtained 
directly  as 


Cf  /Re/ 2  =  a 


(64) 


and 


w  C  /Re 

St  /Ri~  -  - P  -pp  x.  v 

x  (1  -  a) 


(a  -  gc> 


(65) 


NUMERICAL  COMPUTATIONS 


Equations  (58)  and  (59)  are  first  order  differential  equations 
and  can  be  solved  by  any  well  known  integration  technique.  Both 
equations  involve  the  unknowns  "a"  and  a  and  hence  have  to  be  solved 
simultaneous ly . 

The  continuous  blowing  at  the  wall  gives  rise  to  continuous 
boundary  conditions.  These  equations  therefore  are  very  easily 
solved  on  CSMP  (Continuous  Systems  Modeling  Program) .  Numerical 
answers  are  obtained  using  fourth  order  Runge-Kutta,  fixed-step 
integration  technique.  The  step  size  used  was  =  0.0001.  To 
test  the  stability,  this  solution  was  compared  with  the  answers 
obtained  using  rectangular  Simpson's  and  Adam's  method  of  integra¬ 
tion  (all  the  above  techniques  are  available  on  CSMP) .  The  answers 
obtained  by  all  these  techniques  were  within  0.5%. 

The  starting  values  of  "a",  i.e.,  a^n,  are  those  given  in 
Table  I.  For  the  flow  near  the  stagnation  point,  the  "in"  value 
of  "a"  could  be  that  given  in  Table  I  or  II  depending  on  the  mode 
of  injection.  Same  is  true  for  the  initial  value  of  a .  For  wedge 
flows,  the  solution  is  carried  out  up  to  x  =  0.9  ft.-  or  until 
the  value  of  "a"  dropped  below  0.02.  For  the  4:1  elliptic  cylinder, 
the  solution  is  obtained  only  over  the  front  half  while  for  the 
experimental  airfoil,  the  solution  covers  the  entire  surface, 
suction  and  pressure.  Any  time  the  solution  gave  the  value  of  "a" 
below  0.02,  the  computations  were  stopped.  At  this  low  value  of 
the  shear  stress  parameter  the  boundary  layer  is  very  thick  and 
unstable  and  the  prandtl  boundary  layer  assumptions  are  no  longer 
valid. 


The  numerical  calculations  were  carried  out  for  the  following 
values  of  the  variables : 
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RESULTS  AND  DISCUSSION 


To  suppliraent  the  present  analysis ,  the  answers  obtained  from 
Equations  (58)  and  (59)  are  compared  with  the  numerical  computer 
solution  of  Equations  (19)  and  (20).  This  program  was  developed 
by  Davis  and  Werle^2).  The  program  is  extremely  versatile  and 
gives  solutions  for: 

1.  Any  Mach  number,  i.e.,  subsonic  or  supersonic; 

2.  Any  Prandtl  number;  and 

3.  Any  pressure  gradient. 

The  program  uses  a  much  more  realistic  viscosity  -  temperature 
relationship  given  by  Sutherland  formula,  than  the  linear  relation¬ 
ship  used  in  the  analytical  solution.  It  can  be  shown  (see  "Proposals 
for  Further  Research")?  that  if  the  more  realistic  p  a  T0,^5  relation¬ 
ship  between  viscosity  and  temperature  is  used,  the  solution  of  the 
momentum  and  energy  equations  would  be  slightly  different. 

Since  "a"  is  a  strong  function  of  viscosity,  the  differences 
between  the  values  of  "a"  obtained  from  the  analytical  solution  and 
those  obtained  from  the  numerical  solution  (computer  program)  are 
likely  to  be  more  severe  than  the  differences  between  the  two  values 
of  the  wall  temperature.  Also  in  the  computation  of  Stanton  number, 
due  to  the  presence  of  the  factor  (1  -  a)  in  the  denominator,  the 
differences  between  the  two  values  of  Stanton  number  are  magnified. 


Flow  Over  Wedges 


The  effects  of  uniform  blowing  on  the  wall  temperature, 
local  heat  transfer  (in  the  form  of  Stanton  number)  and  local 
shear  stress  parameter  for  flow  over  wedges  are  shown  in  Figures 
2  thorugh  7  (Figure  8  shows  the  velocity  distribution  over  these 
wedges).  Solutions  are  presented  for  ten-fold  variation  in  the 
injection  rate  from  w  =  0.005  to  w  =  0.05.  As  to  be  expected, 
as  the  coolant  injection  rate  increases,  the  wall  temperature 
and  the  shear  stress  at  the  wall  both  decrease.  With  i  .creasing 
distance  from  the  leading  edge,  the  wall  temperature  approaches 
coolant  temperature  (Tc)  and  the  sheai  stress  approaches  zero 
(blown-off  boundary  layer).  At  the  leading  edge,  the  temperature 
gradient  is  infinite  (see  Figures  2  and  5).  This  is  due  to  the 
fact  that  in  the  present  analysis  the  thermal  conductivity  of  the 
wall  has  been  neglected. 

Comparison  between  the  analytical  values  and  numberical 
computer  values  of  the  wall  temperatures  is  very  good.  For  shear 
stress  parameter,  the  numerical  solution  predicts  a  much  steeper 
gradient  than  the  analytical  solution  (for  small  values  of  x)  and 
a  more  gentle  one  for  large  values  of  x.  Agreement  improves  as  the 
injection  rate  increases. 


Flow  at  the  Stagnation  Point  and  Over  the  Surface  of  a  Two-Dimensional 
Body  'in  Crossflow  ’  "  “ 

Effect  of  coolant  injection  through  a  porous  wall  of  4:1 
(ratio  of  major  axis  to  minor  axis)  elliptic  cylinder  is  shown 
in  Figures  9,  10  and  11.  The  velocity  distribution  and  the  pres¬ 
sure  gradient  parameter  X  over  this  cylinder  are  given  in  Figure  12. 
This  figure  also  shows  the  contour  of  the  elliptic  cylinder.  The 
geometric  parameters  of  the  cylinder  are 


.Minor  axis -.2 
*•  Major  axis ' 


(G6) 


U 

e 


(1  +  A)  sin 


(67) 


where 


Minor  axis 
Major  axis 


(68) 


and  6  is  the  angle  measured  in  radians  from  the  front  stagnation 
point. 


x 


0 

/ 

O 


e  cos  (0)  d0 


(69) 


(13) 

These  values  were  obtained  from  Evans  (the  reference> 
length  is  (Major  axis)/2). 

Agreement  between  the  numerical  and  analytical  values  is  very 
good.  Analytical  solution  is  presented  with  and  without  Euler 
transformation.  As  far  as  the  wall  temperature  distribution  is 
concerned,  (Figure  9) ,  the  Euler  transformation  does  not  make  much 
difference.  However,  for  the  shear  stress  parameter,  the  solution 
obtained  with  the  Euler  transformation  agrees  much  better  with  the 
computer  answers  than  that  obtained  without  Euler  transformation. 


The  agreement  between  the  analytical  and  numerical  solutions 
clearly  support  the  validity  of  the  assumptions  made  earlier.  It 
was  assumed  that  for  small  absolute  values  of  N(|n|  <  0.7),  terms 
involving  higher  powers  of  N  (e.g.,  ( e,  dK/d£)2  etc.)  could  be 

neglected  fo>*  the  coefficients  of  n^,  ,  p”7  and  n®.  (In  fact, 

for  n”7  and  n8  all  terms  involving  N  are  negligible.)  Ve  .y  good 
agreement  between  the  numerical  and  analytical  values  shows  that 
the  contribution  by  these  terms  is  indeed  small  enough  to  be 
neglected. 


It  is  necessary  to  realize  that  this  assumption  is  valid  only 
(i)  for  small  injection  flow  rates  (which  are  likely  to  be  en¬ 
countered  in  the  turbine  blade  cooling, ;  and  (ii)  near  the  front 
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stagnation  point. 

Near  the  stagnation  point,  A  is  of  the  order  of  -1.0  and  N 
is  at  least  one  order  of  magnitude  lesser.  Consequently,  the  terms 
involving  N2,  [C(N)^32  etc.  are  very  small.  Away  from  the  stagna¬ 
tion  point  (i.e.,  with  increasing  x) ,  the  absolute  value  of  A 
decreases  while  that  of  N  increases.  Downstream  from  the  crossover 
point  (this  is  the  value  of  x  when  j A J  =  j  N  | )  ,  the  contribution  of 
higher  order  N-terms  becomes  important.  Fortunately  for  most  any 
two-dimensional  body  in  crossflow  for  small  injection  rates,  this 
crossover  point  will  occur  far  enough  downstream  from  the  stagnation 
point  such  that  the  contribution  by  the  n5  or  higher  terms  (of  the 
series)  will  be  small  enough  to  be  neglected.  Hence,  the  validity 
of  the  assumptions  will  hold  for  all  values  of  x. 

As  mentioned  in  the  chapter  on  the  stagnation  point  solution, 
the  initial  conditions  for  the  stagnation  flow  depend  on  the  mode 
of  injection.  If  the  injection  starts  for  x  >  0,  the  initial 
values  are: 


a.  =  1.234  (with  Euler  transformation) 

in 


If,  however,  the  injection  is  assumed  to  initiate  from  x  =  0, 
the  initial  conditions  are  the  values  given  in  Table  II.  These 
values  are  also  presented  in  Figures  13  and  14.  Figure  13  also 
shows  the  stagnation  point  temperature  obtained  by  other  investi¬ 
gators.  The  present  analysis  predicts  slightly  lower  values  of 
the  stagnation  point  temperature  throughout  the  range  of  injection 
considered  (w  =  0.001  to  2  =  0.01).  The  solution  is  valid  for  ary 
two-dimensional  body  in  crossflov;  provided  the  value  of  N  at  the 
stagnation  point  does  not  exceed  (in  die  negative  sense)  -1.0. 

Solutions  of  the  wall  temperature  and  the  shear  stress  for 
the  4:1  elliptic  cylinder  with  the  initial  values  obtained  from 
Table  II  (i.e.,  stagnation  values)  are  also  presented  in  Fig¬ 
ures  9  and  11.  In  all  cases,  the  solutions  quickly  merge  with 
the  one  obtained  with  the  initial  values  from  Table  I. 


Flow  Over  an  Airfoil 

The  stagnation  point  analysis  is  also  applied  to  the  flow  over 
an  airfoil  (snown  in  Figure  15).  The  experimental  velocity  distri-7... 
bution  over  the  airfoil  was  obtained  by  M.  Hussein  and  W.  Tabakoff' 
Readings  were  taken  at  the  room  temperature,  with  UQ  =  142.6  ft/sec 
and  Reynolds  number  (based  on  chord  length)  =  625,000.  The  values 
of  a  and  dA/d£  were  calculated  from  Figure  16  (which  gives  the  velocity 
distribution  over  the  suction  and  the  pressure  side  of  the  airfoil) . 
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present  theory  does  not  give  satisfactory  answer  as  to  whether  the 
iar.inar  boundary  layer  will  undergo  transition  to  turbulent  boundary 
layer  or  will  get  blown  off  when  shear  drops  to  very  low  values.  A 
more  thorough  investigation  is  necessary  to  evaluate  the  effect  of 
injection  on  laminar  boundary  layer  stability. 


CONCLUSIONS 


An  analytical  solution  for  compressible  laminar  boundary  layer 
with  continuous  mass  injection  was  derived.  Here  the  previous  work^1' 
is  extended  such  that  the  solution  is  valid  for  any  pressure  gradient 
including  the  stagnation  point;  for  subsonic  flows  with  a  Pranatl 
number  of  one  and  for  a  linear  viscosity-temperature  relationship. 

The  method,  however,  can  be  broadened  ro  consider  high  Mach  numbers 
and  real  representation  of  Prandtl  number  and  viscosity  temperature 
variation.  Numerical  examples  for  the  wall  temperature  and  shear 
stress  were  obtained  for  continuous  injection  and  the  values  compared 
with  exact  computer  solutions.  Very  good  agreement  between  the  two 
values  of  wall  temperature  was  realized.  The  agreement  in  the  shear 
stress  values  is  fair. 

The  following  general  conclusions  are  drawn: 

(i)  For  transpiration  cooling,  the  wall  temperature  falls 
rapidly  with  increasing  injection  mass  flow  rate  and  the  tempera¬ 
ture  gradient  is  infinite  at  the  leading  edge. 

(ii)  Injection  of  the  coolant  into  the  boundary  layer  makes 
the  layer  unstable.  This  is  most  likely  to  cause  early  transition 
from  laminar  to  turbulent  boundary  layer. 

(iii)  When  the  shear  stress  parameter  "a"  drops  below  0.02, 
the  solution  becomes  unstable. 

(iv)  In  cooling  a  turbine  blade  by  transpiration  cooling 
method,  an  injection  rate  of  0.5%  of  the  main  flow  is  sufficient 
to  maintain  uhe  blade  temperature  below  65%  of  the  mainstream 
temperature . 


(v)  Injection  of  mass  into  the  laminar  boundary  layer  causes 
instabilities  in  the  boundary  layer  and  may  result  either  in  transition 
to  turbulent  boundary  layer  or  laminar  blow-off.  Further  investigation 
is  necessary  to  study  these  phenomena. 
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TABLE  IV 

Explanation  Of  Symbols  Used  In  The  Figures 


Symbol 

w 

Remarks 

O 

o 

0.001 

Each  symbol  in  the  figure 

O 

0.002 

represents  a  data  point  obtain- 

□ 

n 

0.005 

ed  from  the  numerical  computer 

O 

0.01 

solution. 

A 

0.02 

0.05 

(Except  in  Integral  solution,  terms  up 

Figs.  9-22)  to  in  the  series. 

(Except  in  Integral  solution,  terms  up 

Figs.  9-22)  to  r|8  in  the  series  (with 

Euler  Transformation) . 


FIGURE  1  PHYSCIAL  fiOui-.L  FOR  FLOW  OVER  A  WEDGE 


PARAMETER  FOR  THE  TRAHSPRIATION  COOLED  V!, 


TRANSPIRATION  COOLED  WALL, 


TRANSPIRATION  COOLED  WALL  TEMPERATURE  DISTRIBU 


X  FT 

RESS  PARAMETER  FOR  THE  TRANSPIRATION  COOLED  WALL,  A.= 


FIGURE  8  VELOCITY  DISTRIBUTION  OVER  WEDGES 


FIGURE  9 


WALL  TEilPERATU 
CYLINDER 


i 


WITHOUT  EULER 
TRANSFORMATION 


INITIAL  CONDITIONS 
FROM  TABLE  III 
(WITH  EULER  TRANS¬ 
FORMATION) 


WITH  EULER 
TRANSFORMATION 
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DISTRIBUTION  FOR  THE  4:1  ELLIPTIC 


37 


TRANSPIRATION  C 


Nomenclature  for  Figure  13 


Symbol 


■  x 


Reference  j 

Remarks 

15 

Cylinder  in  Crossflow; 

Ue  =  2  sin(x) 

5 

Two-Dimensional  Body  in  Crossflow; 

Ue  =  sin(«x/2) 

6 

Turbine  Blade; 

Velocity  Distribution  obtained  from 

a  figure  in  the  Reference  i 

5 

Same  Turbine  Blade; 

Computer  Solution 

0.2 


CD  00  ■  O  CVI 

O  O'  ^  ^ 


'Ul  A-33W1SIG  lVIif!39 Wd 
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SOLUTION  OBTAINED  WITH 
EULER  TRANSFORMATION 


_i _ L_ 

.a  o.e 


x 


TRANSFER  OVER  THE  SUCTION  SIDE  OF 


0.8 


THE  AIRFOIL 


0  0.2  0,4  0.6  0.8  1,0 

X 


FIGURE  20  WALL  TEMPERATURE  DISTRIBUTION  OVER  THE  PRESSURE  SIDE  OF 
THE  AIRFOIL 
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(Stv^e~)  x  1(3 


12. Or 


FIGURE  21  HEAT  TRANSFER  OVER  THE  PRESSURE  SIDE  OF  THE  AIRFOIL 
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